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Abstract. Generalized #CSP problems include Holant problems with pla- 
narity restrictions; polynomial-time algorithms for such problems include match- 
gates and matchcircuits, which are based on Pfaffians. In particular, they use 
gates which are expressible in terms of a vector of sub-Pfaffians of a skew- 
symmetric matrix. We introduce a new type of circuit based instead on de- 
terminants, with seemingly different expressive power. In these determinantal 
circuits, a gate is represented by the vector of all minors of an arbitrary matrix. 
Determinantal circuits permit a different class of gates. Applications of these 
circuits include a new proof of the Chung-Langlands formula for the number 
of rooted spanning forests of a graph and a strategy for simulating quan- 
tum circuits with closed timelike curves. Monoidal category theory provides a 
useful language for discussing such counting problems, turning combinatorial 
restrictions into categorical properties. We introduce the counting problem 
in monoidal categories and count-preserving functors as a way to study FP 
subclasses of problems in settings which are generally #P-hard. Using this 
machinery we show that, surprisingly, determinantal circuits can be simulated 
by Pfaffian circuits at quadratic cost. 

1. Introduction 

Let Vectc be the category of finite-dimensional vector spaces and linear transfor- 
mations over the base field C. A string diagram [8 in Vectc is a tensor ( contraction) 
network, and fixing such a diagram, the problem of computing the morphism rep- 
resented is the tensor contraction problem, which is in general #P-hard (examples 
include weighted counting constraint satisfaction problems [3]). 

We study complex- valued tensor contraction problems in subcategories of Vectc 
by considering them as diagrams in a monoidal category. For a survey of the rich 
diagrammatic languages that can be specified similarly see [18] and the references 
therein. By a circuit we mean a combinatorial counting problem expressed as a 
string diagram in a monoidal subcategory of Vectc (that is, a tensor contraction 
network). Such diagrams generalize weighted constraint satisfaction problems and 
Boolean circuits; for example fanout and swap are not necessarily available. As 
in the Boolean case, such circuits often have more familiar description languages. 
Subcategories of Vectc faithfully represent Boolean[TT] and quantum circuitspQ, 
counting constraint satisfaction problems, and many other problems[5]. 

Suppose we have a a problem C, for example a counting constraint satisfaction 
problem [5] given by an embedded planar graph with variables or weighted con- 
straints at each vertex. We can consider such a problem to be described by the 
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data of a monoidal world (see e.g. [10j . Chapter 12). We consider an interpretation 
|18j map i: £—> Vectc that gives an interpretation of the problem of interest as a 
tensor contraction problem in Vectc- A second, count-preserving functor h from 
a category C in which contraction is in FP and a subcategory S of Vectc serves 
to characterize the problems which can be solved in polynomial time according to 
a particular contraction scheme. Many such schemes (such as holographic algo- 
rithms |19j ) work by exploiting some combinatorial identity or kernel relating an 
exponential sum (corresponding to performing the tensor contraction by a naive 
algorithm) and a polynomial time operation (often the determinant or Pfafhan of a 
matrix) that yields the same result. They can be viewed as a complementary alter- 
native method to geometric complexity theory |17j in the study of which counting 
problems (such as computing a permanent) may be embedded in a determinant 
computation at polynomial cost. 

We formulate a class of circuits based on determinants and show that it is always 
solvable in polynomial time. This existence of such a class was conjectured in |12j . 
We give some applications of these circuits including a new proof of the Chung- 
Langlands formula [5] for the number of rooted spanning forests of a graph. We 
explain the relationship between Pfaffian circuits |15) (and so matchgates) and 
dctcrminantal circuits and prove a functorial relationship between them. We show 
that, surprisingly, every dctcrminantal circuit can be expressed as a Pfafhan circuit 
at quadratic cost. 

This paper is organized as follows. In Section[2j we describe the counting problem 
in monoidal categories and the setting for our results. In Section [3] we define de- 
terminantal circuits and give applications to the Chung-Langlands rooted spanning 
forest theorem (Section |4| and quantum circuits with postselection-based closed 
timelike curves (P-CTC, Section pif. In Section | wc relate dctcrminantal and 
Pfafhan circuits. 



2. Toward a categorical formulation of counting complexity 

Let M. be a (strict) monoidal category [TJ] with monoidal identity 1» and 
such that Sm = Hovum^m, 1m) is a semiring; call this a semiringed category. 
A monoidal word is a collection of morphisms composed and tensored together to 
form a new morphism [lOj . The problem of counting in a semiringed category M. 
is to determine which morphism in fIom_A/i(ljvf , 1,m) is represented by an arbitrary 
monoidal word in M with \m as its domain and codomain. 

Example 2.1. Consider the monoidal category in which the objects are cartesian 
powers of a two-element set B = {T, F). For example 

{T,F} 2 = { (T, T) , (T, F) , (F, T) , (F, F) } 

which we denote by {T®T, T®F, F®T, F(g)F} since the tensor product of objects 
in this category is the cartesian product of the corresponding sets. Morphisms are 
nonnegative integer valued binary relations such as {2 • (T®F, T®T), 3 • (T®F, T® 
F)}. The composition is F o G = {(/3~f) ■ (a, c) : f3 ■ (a, b) e F, 7 • (b, c) e G} where 
a,/3,7Z;»o- 

This is a semiringed category with Sm = N in which the counting problem is #P- 
hard, and the familiar category of counting constraint satisfaction problems can be 
expressed as a symmetric monoidal subcategory. Every constraint satisfaction prob- 
lem can be visualized by a bipartite graph. One independent set is the variables, 



DETERMINANTAL CIRCUITS 



3 



each of which can be assigned the value true or false. Thus we allow the morphisms 
in Hom(l, B® n ) for each n of the form {1 • (1, F®F® ■ ■ ■ F), 1 • (1, TOT® • • • T)}. 
The other independent set is the constraints. So we allow only those morphisms in 
Hom(i?® Tl . 1) whose coefficients are only zero or one. The string diagrams in this 
category will then be precisely the graphs visualizing boolean constraint satisfac- 
tion problems. In the sequel we generally use the set {0, 1} instead of {T, F} to 
conform with common usage. 

A strict monoidal functor F : M. —* M 1 between semiringed categories is count 
preserving if the induced map F : Sm ~* is a injective morphism of semirings. 
Schemes that generalize holographic algorithms |19j seek a count-preserving functor 
from an category with a FP counting problem to a category with an (in general) 
#P-hard counting problem. 

In each type of circuit, we consider two semiringed categories C and S. Let C 
be a problem of interest. We call C the counting category and S is a subcategory 
of Vectp. Then let / : C — > C be a map that gives an interpretation or encoding of 
the problem as a string diagram in C. By this we mean that for every instance of 
a problem I e C, f(l) is a string diagram that solves this instance of the problem. 

The category C may have a non-intuitive encoding of the problem but has the 
advantage that there exists a polynomial-time algorithm to determine which mor- 
phism of Sm = Hom(l, 1) is represented by an arbitrary monoidal word. We also 
have an interpretation g : C — ► S. Then we want a monoidal functor h such that 
the diagram 




h 



commutes and such that the count is preserved by h. S is the subcategory generated 
by the morphisms in the image of either h o f or g. The induced maps on Sc and 
Ss makes Sc a sub-semiring of Ss- The functor h is called sDet, and sPf for 
dctcrminantal and Pfaffian circuits respectively in the sequel. 

Of course, it is important that the construction represented by the functors is 
implementable in polynomial time. Often this is not a concern, because diagrams 
in C and C are effectively identified, and the problem is expressed in the language 
that will be used to perform the contraction. 

3. DETERMINANTAL CIRCUITS 

Suppose X is a n x m matrix of elements of a field k with rows and columns 
labeled by finite disjoint subsets N and M of N = Zj> - For teN, let Vi = k 2 
be spanned by an orthonormal basis (with inner product) 1^0,^,1 and for finite 
N c N write Vn '■= ®%€NVi- Define the function sDet by 

sDet : Mat fc (n, m) -► V$ ® V M = (C 2 *)®" ® (C 2 )® m 
sDet(X)= ^ det(X JJ )|Z><J| 

I^[n],Jcz [m] 

where |f> = (X), jeAr u< )X (i,j), (J\ = <S>ieM V i,x(*-J) and the indicator function x(h I) = 
if i $ I and 1 if i e J. 
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Figure 1 . An example of a determinantal circuit (wires oriented 
clockwise). The four tensors in Vectc, from left to right, are ob- 
tained by applying sDet to each matrix. Letting V = C 2 , they 
lie in (V*)® 2 <g> V® 2 , (V*)® 2 ®V, V* ® V® 2 , and (V*)® 2 ® V"® 2 
respectively. 



This subdeterminant function sDet induces a strong monoidal functor sDet : 
C — > Vectc from a matrix category to a subcategory T> of Vectc- Let C be the 
free monoidal category described as follows. The objects of C are finite ordered 
subsets of N (which may have repeated elements), with monoidal product union. 
The morphisms arc C-valued matrices with rows and columns labeled by subsets 
of N. If M, N are two matrices with the set of row labels of M equal to the set 
of column labels of N, order them and let N o M = NM be the ordinary matrix 
product, with the resulting matrix inheriting the row labels of N and the column 
labels of M. The monoidal product ®c is the direct sum of labeled matrices. 

Let V be the image of C in Vectc- It will be free dagger symmetric traced [5] 
monoidal subcategory of finite-dimensional C-vector spaces generated by the object 
C 2 , endowed with an orthonormal basis, and morphisms sDet(M) for M a labeled 
matrix. Tensor product and composition/contraction are the usual operations. 

For / : A — » A e C, define tr(/) = det(J + /) and define trace in T> in the usual 
way. For a matrix M representing /, tr(M) = SiWijwi] detMj. 

Proposition 3.1. C is a strict dagger symmetric monoidal category. 

Proof. We have already defined ® for C. We need to show that it is a bifunctor. 
For A cz N, id^ is the identity matrix with row and column labels A. It is easy to 
see that for any 4,B c N, id^tguds = id^B- Now for morphisms W,X,Y } Z e 
Mor(C), W ®X oY ® Z = {W ® X){Y ® Z) = WY®XZ = (W o Y) ® (X o Z), 
so ® is indeed a bifunctor. 

For A,B,C e Ob(C) , the associator a A BC ■ {A ® B) ® C -> A ® (B ® C) is just 
equality by the associativity of matrix direct product. The unit for C, denoted 1, 
is the empty set. Then : 1 ® A — > A and pa ■ A ® 1 — * A are also equality since 
it is union with 0. It is clear that a, A, and p are natural isomorphisms. 

We need to check that the diagrams from MacLane's Coherence Theorem com- 
mute. First let us check, for A, B e Ob(C): 

(A ® 1) ® B — A ® (1 ® B) 

id A ®A B 

¥ 

A®B 

(A®t)®B = (A u 0) u B is mapped to A u B by pa ® ids via equality. Then a 
maps (A u 0) u B to A u (0 u B) via equality. This is then mapped to A u i? by 
id^^As via equality, and the diagram commutes. 
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Now let us check the second diagram, for A, B,C,D e Ob(C): 



((C®A)®B)®D- 



(C®(A®B))®D 



(C®A)®(B®D) 



C®(A®(B®D)) 



id c ®ct 



C®((A®B)®D) 



The object ((C®A)®B)®D) = ((Cu A) u B) u D is mapped toCu(Au(BuD)) 
by (idc ®a)o(a)o(a®idi)) via equality. Similarly, it is mapped to Cu (Ad (BuD)) 
by a o a via equality. This diagram also commutes and so C is a monoidal category. 
Furthermore, since a, A, and p are equalities, C is a strict monoidal category. 

The braiding for C is a map ca,b : A® B — » £? (g) A, A, B e Ob(C). It is given by 
the matrix 

£ A 

a ( o r 

CA < S = S (, 1 

We need to check that the following diagrams commute for A,B,C e Ob(C): 



(B®A)®C- 



B®{A®C) 



Jd B ®C A ,C 



{A®B)®C 



B®{C®A) 



A®{B® C) A -^S-\B ®C)®A 



{A®B)®C 



(B®A)®C- 



B®{A® C) 




B® (C® A) 



A®(B®C) ±^\b ®C)®A 




The first diagram commutes by noting that 
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The second diagram commutes since c~j^~ A = ca,b (which implies the category is 
symmetric) for any A, B e Ob(C) and so the second diagram is the same as the 
first. 

The dagger for C is given by matrix transpose and the identity on objects. Clearly 
id^ = id^ = id A - Given X, Y e Mor(C), X : A -> B. Y : B -> C, {X o F)t = 
(XY) T = Y T X T = y t ol':C-»A Lastly JSftt = X TT = X. 

We also need the dagger to satisfy two extra properties since we are working 
in a monoidal category. First, given X, Y e Mor(C), (X © Y") 1 " = (X © Y~) T = 
X T © Y T = X^ © FL Secondly, a, A, and /? should all be unitary (its inverse is 
equal to its dagger). Since they are all the identity morphism, this is also satisfied. 
Thus C is indeed a strict dagger symmetric monoidal category. □ 

Theorem 3.2. The map sDet defines a strict monoidal functor which is an equiv- 
alence (in fact, an isomorphism) of dagger symmetric traced categories. Thus while 
computing a trace in Vectc is in general ^P-hard, in the image of sDet it can be 
computed in polynomial time. 



We prove this in two parts as Lemmata 3.3 and |3.4 



Lemma 3.3. The map sDet defines a strict monoidal functor which is an equiva- 
lence of monoidal categories. 

Proof. We have defined how sDet acts on matrices, and so morphisms in C. It takes 
an object A e Ob(C) to sDet(^4) = V A = ^aVi. 

First we must show that sDet is a functor, i.e. that it respects composition and 
that sDet^d^) = id sDct ( J 4). Suppose X e Home(I, J), Y e Homc(J, K) so X is 
a matrix with row labels /, column labels J and Y has row labels J and column 
labels K: 

sDet(F) o sDet(X) = E E E det(X tf ) det(y ifc )l»X*l 

= ^ E ^{*Y lk )\i){k\ = sDet(Xr) 

where the middle equality is the Cauchy-Binet formula. Now in C, id a is the identity 
matrix with row and column labels A. Then sDct(id^) = YjiaA l-OC'l which is the 
identity morphism for the object sDet(A) in T>, and sDet is indeed a functor. 

For sDet to be a monoidal functor, we must demonstrate two additional prop- 
erties. First we must show that sDet(A © B) = sDet(A) ©sDet(_B). Let / and J 
be the rows and columns of A, respectively. Let I' and J' be likewise for B. A 
straightforward calculation gives 



sDct(A©B)= J] det(A@B) uv \UXV\ 

= ^ E det(A UnIiVn j)det(B VnPiVn j,)\Unr)\UnI'XVn J\(V n J'\ 

= E E E E det(A uv )det(B u , v ,)\U)\U'XV\(V'\ = sDet(A) © sDet(B). 
ye/ c/'q/' vcj ycj' 

Secondly we must show there are morphisms Fo : 1-p — > sDet(lx>) (the unit in T> 
is the base field C) and for any A, B e Ob(C), f\ : sDet(A)®sDet(S) sDet(A©73) 
satisfying certain axioms expressed as commutative diagrams. 
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Since sDet(0) = (g«=0Vi = C, F is simply equality. Similarly for objects A and 

B, 

sDet(A®B) = sDet(A uB) = ^Aui?^, = (g) (®j e s^) = sDet(A) <g> 

sDet(i3); so F\ is equality. In the following diagrams, we shall call sDet simply F. 

Let a', A', p' be the natural transformations for V. Note that all three are equal- 
ities. For A, B,C e Ob(C), the following must commute: 



F(A) <g (F(B) (g F(C)) 



F(A)(g) (F(5<gC)) 
F(A<g> (5(g)(7)) 



F(a) 



(A) ® F(B)) (g) F(C) 

F 1 ®id J?(c) 

(F(A®B)®F(C)) 
-^-f((A®B)®C0 



(g) 1' 



id F (B) ®F 



F(fl)®F(l) 



F(B®1) 



1'0F(B) 



F{B) 



F (x)id F(B) 



F(1)®F(B) 



F(A) 



F{t®B) 



The diagrams trivially commute as all of the maps are identities. So sDet is a 
strong monoidal functor. Since F , F\ are equalities, it is in a fact a strict monoidal 
functor. 

Lastly, we want to say that C and T> are equivalent as monoidal categories. By 
definition of T>, sDet surjects onto objects and morphisms, so it is a full functor. 
Now consider Hom(A, B) for objects A, B e Ob(C). Let X e Hom(A, £?). sDet(X) 
contains all the entries of X as coefficients in the sum since the entries of X are lxl 
minors, and X is determined by its image sDet(JT). Thus sDet induces an injection 
on Hom(j4, B) — » Hom(sDet(A), sDet(B)), and the functor is faithful. Thus it is 
an equivalence. However, it is not quite an isomorphism as sDet does not give a 
bijection on objects. □ 

We have yet to define the braiding and dagger for T> required to state Theorem 
|3.2| For sDet to respect the braiding, we need the following diagram to commute: 



F{A)®F{B) 



Fi 



F(A (g) B) 



F(B)®F(A) 



Fi 



F{ca,b) 



F(B (g) A) 



Since F± is the identity, we just need that CpiA),F(B) = F(ca,b)- F{ca.b) = 
|00><00| + |01><10| + |10><01| - |11><11|. We define the braiding for V to be this. 
For the dagger, consider X e Mor(C) with row labels / and column labels J. 

BDet(Xt)= J] det(X?)K>01 = 
^ det(Jfy)|jX»l = sDet(X) T 



8 



DETERMINANTAL CIRCUITS 



So the dagger for V is the normal dagger in Vect^. 

Again for / : A — > A e C, define tr(/) = det(7 + /) and define trace in V in the 
usual way. 

Lemma 3.4. The map sDet defines a strict monoidal functor which is an equiva- 
lence of dagger symmetric traced categories. 

Proof. By construction, sDet respects the braiding. We also showed that this func- 
tor respects the normal dagger for linear transformations. Theorem |3.6| and Propo- 
sition 3.7 below shows that sDet induces the identity map from Hom(le,le) — * 
Hom(lx), Id) and thus respects the trace. 

□ 

Remark 3.5. This braiding is not the usual braiding for Vectc- Thus while the 
functor sDet is count-preserving, the count will not be the same as if the standard 
braiding is used. 

Using the operations of © and matrix multiplication, we can transform any 
string diagram in C into a diagram with a single matrix, M, and thus evaluate the 
dctcrminantal circuit efficiently. 

Theorem 3.6. The time complexity of computing the trace of a determinantal 
circuit in C is 0{dw u} ) = 0{dw u} + c") where d is the depth of the circuit, w is 
the maximum width, c is width at the input and output (so can be chosen to be the 
minimum width), and u is the exponent of matrix multiplication. 

Proof. We have an n x n matrix with equal row and column labels, which we may 
assume to be 1, . . . , n. Then 

sDet(M)= J] det(M ItJ )\I)(J\ 

/,jQ[n] 

and contracting this against itself gives 

2 det(M IJ )(J\I)(J\I)= £ det Mi,/. 

I,Jc[ n ] fc[n] 

That is, the trace of a matrix M in C is the exponentially large sum of its 2™ 



principal minors; we claim that det(7 + A) is precisely this sum (Proposition 3.7). 
This is a classically known fact but we include a proof here for completeness. This 
enables us to compute this number in time n". □ 

The following identity is well-known (e.g. it can be derived from results in [7]); 
we include a proof for completeness. 

Proposition 3.7. Given two n x n matrices TV, M, with N diagonal, 

det (AT + M) = Yj det(Mi) det(7V / ) 
iqM 

where Mj = Mj j and Nj is the minor of N formed by removing the rows and 
columns labeled I. In particular, det(7 + M) = 2ic=[n] det(Mi). 

Proof. Let Ui be the columns of M and Vi the columns of TV, i = 1, . . . , n. Then 
det(TV + M) = A"=i ( u i + v i)- N° w let /3i, . . . , /3„ be variables and 7 : /3j >-» {v i7 ttj} 
be an assignment of either Vi or ttj to every variable . There are of course 2 n such 



DETERMINANTAL CIRCUITS 



9 



assignments. Summing over all possible 7, we can expand the above into the sum 
of "simple wedges" : 

n 

det(iV + A/)=2/\ 7 (ft) 

7 i=l 

Let us fix a particular 7. Let W be the matrix whose i th column is 7 (/3i). Then 

n 

det(W) =/\ 7 (A) 
i=i 

Let J7 = {i : ft = Ui} ordered in increasing order. Let ?7 = {i : Pi = v^}, also 
ordered in increasing order. Then Wjj = Mjj and Wfj = N-q. 

Now we must show that det(W / ) = det(Mjj) det(Njj). Since U is an ordered set, 
there is no ambiguity in letting U(j) be the j entry of U. Consider the matrix 
transformation where the element Wutj\u(j) l—> Wjj. This is done by swapping 
rows j and U(J) and then swapping columns j and U{j). Since we do two swaps, 
this transformation does not change the determinant; do this for all j e [\U\]. Then 
perform the matrix transformation Wf}rj\ xjtj\ ' * Wjj for all j e [\U\]. This yields 
the following matrix: 

~M V " 

X Nfj_ 

Here X is the portion of columns of M specified by U that don't lie in the minor Mjj 
(under the above permutation). Since N is diagonal, the entries of the columns of N 
specified by U not in the minor Nfj will always be 0. The determinant of this matrix 
is det (Mu) det(ATp),. Summing over all of assignments 7 yields the result. □ 

Note that while V could be equipped with the object duality structure (A, A* , ia^a) 
from the category of finite-dimensional vector spaces to obtain a dagger closed com- 
pact category, the matrix category C is not a closed compact category: it lacks the 
morphisms %a and ca- The morphism ca ■ A (x) A* — ► I would have to be the sDet 
of a 2 x matrix, or the composition of several morphisms to obtain one of this 
type. 

Proposition 3.8. C does not have duals for objects. 

Proof. We cannot have = sDet(M) for any M. The morphism we want is 
|00> + |11>, but there is a unique 2x0 matrix M and sDet(M) = |00>. 

□ 

As a consequence, we really do have to work with traced categories [H] rather 
than the more convenient dagger closed compact categories. 

A diagram in the equivalent categories C,T> is called a determinantal circuit] 
when the morphism represented is a field element, it computes the partition func- 
tion, i.e. counts the weighted number of solutions to the #C5P it represents. 
Because these categories have a traced, dagger braided monoidal category struc- 
ture, they come with a corresponding graphical language [TS] . 

It is also a question of interest which tensors are determinantal. One can test 
whether a vector can be the set of determinants of minors from a matrix using the 
Pluckcr relations to obtain the relations among general minors of matrices. On the 
other hand, for minors of a fixed size this is an open problem [2]. 
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4. Applications 



4.1. Multipaths. We now discuss a specific diagrammatic language, and a way 
of viewing what it is a determinantal circuit counts in terms of multipaths. This 
point of view makes it easier to apply determinantal circuits to specific counting 
problems. 

Suppose we have a determinantal circuit. Our convention shall be that objects 
will be represented by arrows moving from left to right. Vertical stacking will de- 
note tensor product, ordered from top to bottom. In our circuit, we can group the 
morphisms into a sequence of stacks, where each stack is simply a group of mor- 
phisms which are aligned vertically. Each stack is also a single morphism obtained 
by tensoring the morphisms in the stack. 

Consider the k th stack in a determinantal circuit and suppose the incoming wires 
are labeled with the set Ik- For c J t) let \i k be the binary string bi k 1 bi k 2 • • • bj k 
where Ik. m is the m th element of Ik (ordered from top to bottom in the diagram), 
y = and where bi k = 1 if Ik,m 6 an d otherwise. The outgoing wires are 
labeled with the set Ik+i and for ik+i ^ Ik+i- The k th stack can be written as a 
morphism in T> as 



where Wi ht % k+1 is the appropriate scalar. 

Now, for a fixed i k c / fe) i k+1 c J fe+1) we say that W ikiik+1 \xi k XXi k+ i I is a 
(weighted) multipath passing through the k th stack. It specifies which incoming 
wires and outgoing wires of the k th stack are activated. Note that unless \ik\ = 
|*fc+i|, Wi kt i k+1 = as it will be the determinant of a non-square matrix, justifying 
the term "multipath". 

Now we say that a multipath through a determinantal circuit with n stacks is a 
tensor of the form 



for fixed sets i\, . . . , i n . If the determinantal circuit is collapsed to a matrix M, 
sDct(M) is the sum of tensors of the form 



However, unless i 2 = i' 2 ,--.,i n = i' n , this tensor will be zero. Also, given the 
cyclical invariance of the trace, it must be that i\ = i\ as well. Thus the tensor 
represented by sDet(M) is the sum of weighted multipaths of the circuit. Taking 
the trace then sums up the weights. 

Thus the value of a determinantal circuit is the sum of the weights of its mul- 
tipaths. The definition we have given for multipath translates into a collection of 
paths, in the traditional sense, in the given determinantal circuit considered as a 
graph. 

4.2. Recovering the matrix tree theorem and variants. We describe how to 
recover the Chung-Langlands rooted spanning forest theorem [5] by expressing it 
as a determinantal circuit. The theorem tells us that if B is the incidence matrix 
of a graph G with an arbitrary orientation placed on it, then the number of rooted 
spanning forests is det(I + BB T ). An example of a graph is given in Figure]^ a) 
and the determinantal circuit constructed for it in Figure [2](b) . Our proof differs 
from that of Chung and Langlands in that we give an explicit interpretation of 
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every permutation in the calculation of the determinant of a principal minor as a 
structure in the graph. 

Choose an arbitrary orientation on G, and consider a string diagram in C which 
is equipped with a C-morphism (node) for every edge and vertex of G. A node 
corresponding to an edge is connected to a node corresponding to a vertex if the 
edge and vertex are incident in G. This edge is oriented toward the vertex node 
if that vertex is a sink for the edge in G; otherwise the edge is oriented towards 
the edge node. Note that this orientation is not the orientation reflecting primal 
and dual objects in the category. It is simply a tool to help explain the correct 
interpretation of the circuit. Let us call this portion of the diagram Z. It is portion 
of the circuit inside the dashed box in Figure [2jb). 

To express the problem as a determinantal circuit, we will associate matrices 
(morphisms of C) with the nodes of Z . Then applying the simplification as in the 
proof of Theorem |3.6| using the operations of © and matrix multiplication, we will 
obtain the matrix B. Let us imagine that the edge nodes of Z are arranged in 
a vertical fashion. To the right of them, we have the vertex nodes also arranged 
vertically. From top to bottom, vertex nodes will take values 1,2, ... . Similarly, 
edges nodes will take values a, b, . . . . 

Now in Z, for every edge node, let there be a single edge extending to the left. 
Label this edge a if it comes from a, b if it comes from b, and so on. Next we label 
the edges from edge nodes to vertex nodes. Let the edges incident to 1 be labeled 
li, I2, • • ■ , l n where n is the degree of 1 in G. Do likewise for the node 2 and so 
on. Finally extend a single edge from every vertex node rightward. 

For every edge node e, associate a 1 x 2 matrix with it, M e . Let the entry in the 
first row and column of M t be —1 if the first column label denotes an output edge 
oriented away from e. Otherwise, let it be 1. Define the other entry of M t similarly. 
So M t = [1 - 1] or [-1 1]. 

For every vertex node v, associate a n x 1 matrix with it, M v . The row labels 
in order from top to bottom are V\, v-x, . . . , v n where n is the degree of v in G. Let 
the output label be the number v. Let M v have all entries 1. 

Although in general wc suppress it in pictures, whenever two wires cross, we put 
the braiding matrix (as discussed above) on the crossing. 

Lemma 4.1. Using the operations of matrix multiplication and ©, the matrices in 
Z collapse to the incidence matrix of G with some orientation placed on it. 

Proof. Let E be the matrix formed by taking the direct sum of all the matrices 
on the edge nodes. Let V be the direct sum of all the matrices on vertex nodes. 
A = EV will be the matrix that the morphisms in Z collapse to. 

The rows of E are labeled with the edges of G and the columns of V are labeled 
with the vertices of G. Let e be an edge of G and v a vertex. Let A ev be the 
element of A in the e th row and v th column. 

The e th row of E contains a +1 in the column Vi and +1 in column Vj (where i 
and j are two distinct numbers between 1 and the degree of v). If the v th column 
in V does not contain a in both the v\ h and u*' 1 row (the case where e is not 
incident to v), then it will contain a 1 in either the v\ row or the v* h row, but not 
both as there are no two nodes in Z that have two edges in between them. In that 
case v is incident to e in G and A ev is +1. 

So A ev is +1 if and only if e is incident to v in G and otherwise is 0. Now suppose 
that A ev = 1. Let v' be the other endpoint of e. We need to say that A ev > = — 1. 
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Now let Vi be the edge in the diagram that connects the node corresponding to e 
to the node corresponding to v. Then the e th row of E has a 1 in the vf 1 column. 
Then let v'a be the edge in the diagram that connects the node corresponding to e 
to the node corresponding ton v' . Since the e th row of E has exactly two nonzero 
entries, one being 1 and the other being -1, the v'jth column must have a -1 in the 
e th row meaning that A ev > = — 1. 

So indeed A is an incidence matrix of G with some orientation placed on it. 

□ 

Now we reflect Z across a vertical line and take the transpose of all matrices 
associated with nodes of Z. We call this Z T . Indeed, Z T collapses to the matrix 
B T . The outputs of Z are precisely the inputs of Z T and vice versa. So we glue 
Z and Z T together to get the circuit we want. In Figures [2^ a) and (b), we give an 
example of a graph and its transformation into a circuit. Figure [2^c) is equivalent 
to[2|h) but we have multiplied our two stacks of vertex nodes. 

Let us turn our attention to the kinds of multipaths our circuit admits. Consider 
sDet(M e ), where M € is the matrix in Z associated with the edge node e. It is of 
the form |0)(00| + |1)(10| + |1)(01|. This says that for any multipath, e cannot have 
two activated output edges. However, it may not have any path passing through it 
at all. The signs of the various tensors denote the weight of that part of the path. 

Let M v be the matrix associated with the vertex node v (vertex nodes as in the 
collapsed version in Figure [2jc)). Let n be the degree of the corresponding vertex 
v in G. Since M v is a rank one matrix, the only nonzero principal minors are lxl 
minors and the empty matrix. So 

n 

sDet{M v ) = |0 Vl • • • 0^X0^ • • • Vn | + £ |x 4 >< Xl | 

»=1 

where \xi) = {0^ • • ■ \ Vi ■ ■ ■ 0u„) and (xi\ defined similarly. This means that for any 
valid multipath, a vertex node either has exactly one activated input edge and one 
activated output edge or no path through it at all. It is not possible for two are 
more separate paths in the multipath to pass through v. 

Lemma 4.2. Every multipath of ZZ T is a subgraph ofG. 

Proof. Here we are considering multipaths without their weights to show that every 
multipath of ZZ T may be associated to a subgraph of G. For a particular multipath 
we say that an edge node in Z is activated if there is a path in the multipath that 
passes through it. Similarly for edge nodes in Z T . For this multipath, let n(E)z 
be the set of activated edge nodes in Z and (i(E) z t the activated edge nodes in 
Z T . It must be that [i(E)z = li{E)z T ■ Then a multipath tells us which edges arc 
activated in our subgraph, and we include those vertices incident to each of these 
edges. 

□ 

In general, there may be more than one multipath that is associated with a par- 
ticular subgraph of G. Furthermore, these different multipaths may have different 
weights. 

First, let us consider the weight of a path in our multipath. If we consider our 
circuit to have three stacks (as in Figure [2jc)), let us consider the part of the path 
from the first stack to the second, pi . If the subpath travels along an edge oriented 
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towards an edge node, the subpath has weight +1. If the edge is oriented towards 
a vertex node, the subpath has weight — 1. This follows from the way we defined 
the matrix M e for an edge node e. The subpath from the second stack to the third 
stack, p 2 , satisfies the same properties. 

Then the weight of the path is the weight of p\ times the weight of p 2 . So if 
a path has its two subpaths either both oriented towards a vertex node or both 
towards edge nodes, the path has weight +1. If one of the subpaths is oriented 
towards a vertex node and the other subpath is oriented towards an edge node, the 
weight of the path is —1. Now we say that the unassociated weight of a multipath 
is the product of the weights of its constituent paths. 

Recall that a multipath \i must have \i(E)z = t l (E)z T ^ so we can view \i as a 
permutation of the edges in G associated with the edge nodes in p{E) z . This is 
well defined as two different paths cannot emanate from the same edge node in Z 
and thus two different paths cannot both pass through the same edge node in Z T . 

Define the weight of a multipath to be its unassociated weight times the sign 
of this permutation. The reason for this is that every time two wires cross in the 
multipath, we put the braiding matrix at the crossing. This will flip the sign of the 
unassociated weight. 

Definition 4.3. A straightline multipath is one were every path passes through 
ez and e z r for some edge e. By e z , we mean the edge node associated with e in 
Z. e Z T is the edge node associated with e in Z T . 

Lemma 4.4. Every straightline multipath has weight 1. 

Proof. Let fi be a straightline multipath. Viewing fj, as a permutation, it is obviously 
the identity permutation so the weight of /x is the unassociated weight of p. Now 
let us consider an arbitrary path, p, in /x. Let p connect the edge nodes ez and e z r 
by passing through vertex node v. Since ZZ T is symmetric across a vertical line, 
the subpath from ez to v and the subpath from v to e Z T are either both oriented 
towards v or both oriented towards their respective edge nodes. So this path has 
weight 1. Since the path we chose was arbitrary multiplying all the constituent 
paths together gives that the weight of the multipath is 1. 

□ 

Lemma 4.5. Every spanning forest is represented only by straightline multipaths. 
Furthermore every spanning forest has a representation as a multipath. 

Proof. Since any spanning forest can be broken into the disjoint union of trees, it 
is sufficient to show trees only have straightline multipath representations. So let 
us suppose that there is a multipath /j, that is not straightline that represents some 
tree. Since /i is not the identity permutation on the edges activated, we can write 
it as the product of disjoint cycles. Let us consider one of these cycles. This will 
be represented by a submultipath fi' with m paths. Then since fjf is cycle, without 
loss of generality we can represent the paths Pi in // as Pi = e(l) z — » e(2) Z T, 
P 2 = e(2) z — » e(3) Z T, ■ ■ ■ ,P m = e(m) z — > e(l) z r However, this means that e(l) 
and e(2) share a vertex, e(2) and e(3) share a vertex, and so on. So [il represents 
a cycle in G, and /i cannot represent a tree. 

The second part of the lemma is obvious. In any spanning forest, there is a 
way to associate any edge with a vertex adjacent to it such that no two edges are 
associated with the same vertex. Then consider a straightline multipath where each 
path through an edge node travels through the associated vertex node. 
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□ 

Lemma 4.6. If /i is a multipath which is a cyclic permutation, its weight is -1 

Proof. In the previous lemma, we showed that \i represents a cycle in G. Now 
suppose that the orientation we put on G was such that the arrows travel around 
the cycle. Let e and e' be two edges in this cycle that share a vertex v. In fi, there 
is either a path from ez to e' ZT or e' z to e Z T. Suppose the path is from ez to e' ZT . 
The subpath from ez to v in ZZ T has the opposite weight of the subpath from v 
to e' ZT . So the path between the two has weight -1. If the path had been from e' z 
to ezT, this would still be true. 

Now suppose we flip an arrow on edge e. Then let e~ be one edge sharing a 
vertex with e and e + the other. Then the weight of the path from e z to e Z T (and 
from ez to e^ T ) would flip sign. So would the weight of the path from ez to e J T . 
So the unassociated weight of [i will remain the same. Since we can get to any 
orientation on this cycle by flipping arrows, the unassociated weight of this cycle is 
orientation independent. 

Every path in n has weight -1. So the unassociated weight of fi is (— l) k , where 
k is the number of edges of the cycle. Now if k is odd, fi has unassociated weight 
-1. The sign of the permutation of /i is the sign of permuting an odd number of 
elements cyclically. The sign of such a permutation is 1. So the weight of \i is -1. If 
k is even, the unassociated weight of /i is 1. The sign of permuting an even number 
of elements cyclically is -1. So /i has weight -1 once again. 

□ 

Thus we obtain our version of Theorem 2 of [5] as follows. 

Theorem 4.7. The value tr(ZZ T ) = det(I + B T B) is the number of rooted span- 
ning forests of G. Furthermore, det(/a; + BB T ) gives a polynomial where the coef- 
ficient of x k is the number of rooted spanning forests of G with k roots. 

Proof. First let us consider a straightline multipath that represents a spanning 
forest in G. We look at the submultipaths representing the disjoint trees. If a 
multipath fi represents a tree with t edges, there are t paths in the fi. So one of the 
t + 1 vertex nodes associated with a vertex in the tree does not have a path going 
through it. This will be the root of the tree. 

Since fi is straightline, it has weight 1. So every multipath representing an rooted 
spanning forest contributes a +1 to the sum. It is not hard to see that there is 
exactly one straightline multipath for every rooted spanning forest and this is the 
only way to represent a particular rooted spanning forest making use of a previous 
lemma. We now need to argue that there are multiple ways to represent subgraphs 
with cycles and that their weights will cancel out. 

It turns out that every cycle has four multipath representations and that two have 
weight 1 and two have weight -1. This will be sufficient for proving the result. Let 
us consider a cycle in G with k edges {e(i)}. Then we can represent this cycle with 
paths Pi : e(l) z - e(2) Z T,P 2 : e(2) z -+ e(3)#r, . . . , P k ■ e(k) z -+ e(l) Z T which 
will have weight -1 as a multipath by the previous lemma. We can also represent it 
with P{ : e(2) z - e(\) Z T,P' 2 : e(3) z - e(2) Z T , . . . , P' k : e(l) z - e(k) Z T which is 
the previous multipath reflected across a vertical axis. This will also have weight 
-1. 
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(c) Equivalent circuit 



Figure 2. Transforming a rooted graph to a determinantal circuit. 



These are the only two ways to represent a cycle in G where the permutation 
induced by the multipath is not the identity. This is because a nontrivial per- 
mutation tells us how to walk around a given subgraph and there are only two 
ways to walk around a cycle: clockwise and counter-clockwise. However, there will 
be two representations of a cycle with straightline multipaths. Consider what a 
straightline multipath is. It tells us which edges are activated and associates to 
each edge a vertex incident to it. In a cycle, there is only two ways to make such 
an identification. 

Since we have two straightline multipaths representing the cycle, +2 is con- 
tributed to the overall value for every cycle. But of course this is canceled out the 
by the -2 contributed by the other two multipaths representing the cycle. Now if 
a multipath contains a submultipath that is a cycle, taking that submultipath and 
replacing it with one of opposite sign gives a multipath of the same subgraph with 
a different sign. So the only multipaths whose weights actually contribute to the 
sum are those representing rooted spanning forests. 

In the above setup, the Laplacian would actually be Z T Z. However, note that 
in the circuit, we can move the Z T portion around so that the circuit has Z T to the 
left of Z. This is to say that the trace of the category satisfies the trace property: 
det(7 + XY) = dct(7 + YX). The multipaths will be the same (albeit deformed) 
and will still retain the combinatorial meaning discussed above. When considering 
a minor of BB T of size k, that corresponds to a multipath in Z T Z with k paths. 
We can think of this as activating k vertices in G. 
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We have already shown that det(I + B T B) = det(I + BB T ) is the total number 
of rooted spanning forests. Now consider 

det(Ja; + BB T ) = ( 2 BBf\x k 

fe=0 ^ JQ[|V|] ' 
\I\ = \V\-k 

So a minor of BB T of size \V\ — k is one where exactly k of the vertices don't have 
a path passing through them. By our correspondence, this means those k vertices 
are roots in the rooted spanning forest. So indeed this polynomial is a generating 
function for the number of rooted spanning forests where the exponent denotes the 
number of roots. 

□ 

Now we discuss how to prove Kirchoff's Matrix Tree Theorem as a corollary. 
However, this is the statement about the cofactors of the Laplacian of a graph. 
So we once again consider the determinant al circuit Z T Z which gives the trace 
det(I + BB T ). 

Suppose we have a determinantal circuit that collapses to a matrix M. Then take 
the determinant of some principal minor M' . If M' has labels /, this corresponds to 
the value of activating just the external edges (those that wrap around the circuit) 
labeled I in the circuit. Then taking det(J + M) is just summing up the values 
over all possible subsets of the external edges. Since the matrix tree theorem deals 
with cofactors, in our circuit Z T Z, we will consider multipaths with all edges in 
the diagram turned on except for one. 

Corollary 4.8. (Kirchoff's Matrix Tree Theorem) The number of spanning trees 
of a graph is the absolute value of any cof actor of its Laplacian. 

Proof. Consider the absolute value of an arbitrary cofactor of Z T Z. This will be 
the value of Z T Z will all the edges activated except for one. This corresponds to 
multipaths where all vertex nodes of our graph are activated except for one. Let 
the activated vertices be denoted by V. We only need to consider straightline 
multipaths as all other multipaths will be canceled out. 

Take a given straightline multipath. For every vertex node, there will be an 
associated edge node. So if \V\ is the number of vertices in our graph, we have an 
rooted spanning forest with \V\ — 1 edges. This is a spanning tree. So the number 
of rooted spanning forests on V' is the number of spanning trees of our graph. 

What is happening here is that we are specifying a root for our spanning trees 
(the vertex omitted). Which vertex we specify as the root doesn't matter. Then 
we are counting maximal spanning forests with that root. These will simply be 
spanning trees. 

□ 

4.3. Simulating quantum circuits in the presence of closed timelike curves. 

Determinantal circuits define a class of tensor networks with a polynomial-time con- 
traction algorithm. An immediate consequence is that certain types of quantum 
circuits (or more generally tensor networks possibly including preparations and 
postselection) can be simulated efficiently using this technique. Essentially these 
are the tensor networks of the type shown in Figure [I] (with arbitrarily many wires 
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and transformations): a block of multilinear logic with the same number of inputs 
and outputs, with those inputs and outputs connected via nested loops. 

The loop in such a circuit corresponds to a postselected closed timelike curve 
(P-CTC) [T3]. These are constructed from bell states (cups) and postselected bell 
costates (measurements that act as caps) . The resulting logical category of circuits 
represent physical experiments (which, if they contain an embedded contradiction, 
have count zero |16]). 

If postselection and preparation operations beyond those necessary to construct 
the cap and cup are not allowed, all the transformations represented by the boxes 
are restricted to be unitary. If sufficient additional postselection or preparation 
is allowed, arbitrary transformations of the form sDet(>l) can be appear in the 
boxes. The resulting linear transformation can be interpreted physically in terms 
of a "quantum SVD" as in |16j , where the amount of postselection or preparation 
required to represent an arbitrary linear transformation is also discussed. 

5. Relation to Pfaffian circuits 

Pfaffian Circuits l'ij were introduced as a reformulation of matchcircuits [TU] . 
We present an slightly different definition using category theory. We want to know 
what the relation of determinant circuits is with respect to Pfaffian circuits. The 
most natural way of doing this is functoriality. 

We now define the category that gives us Pfaffian circuits. Consider the set 
J% x {0, 1}, where ^ is the set of labeled skew-symmetric matrices. Furthermore, 
the columns and rows should have the same labels in the same order. The label sets 
are subsets of N. As before, for ieN, let Vi = C 2 be spanned by an orthonormal 
basis (with inner product) t^o, Uj,i and for JV c H write Vm '■= ®isNVi- Now let us 
consider the following function: 

sPf : J( x {0,1} ->V£®V N 

sPf(Af,0) = 2 Pf (Mj) | J> 
;cjv 

sPf (M, 1) = J] Pf (Mj)(I\ 

where |f> = (X) ieJV u^ x (i,j), (J\ = ® ieM v * x (i,j) and tne indicator function x{h -0 = 
if i $ I and 1 if i e I. Here Mj means the principal minor of M with row and 
column labels I. Mj means the principal minor of M with the rows and columns 
labeled I removed. We will use the convention that sPf(M,0) will be denoted 
sPf(M) and sPf(M, 1) will be denoted sPf v (M). 

The sPf function lets us define a monoidal subcategory of Vectc- Let 3 s be 
the free monoidal category defined as follows. The objects are of the form Vjy for 
ordered subsets of N, the tensor product being the usual one. The morphisms of 
are generated by elements from the image of sPf . Composition and tensor product 
will be inherited from Vectc- 

Theorem 5.1. 9* is a strict monoidal category with daggers. 

Proof. By our definition of 0*, it will be the smallest monoidal subcategory of Vectc 
containing the generating morphisms with the specified objects. Furthermore, we 
can assume we are dealing with the strict version of Vectc without loss of generality. 
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So the a, A, and p maps that & inherits will be identities. We want to show that 
the identity morphism is actually generated by our specified morphisms. Consider 
the following matrix for an object A: 



Ia 



A 
A 



A 


A 





1 


-1 






LetL A = sPf(/ A ) = \0 A Q A )+\l A \ A )uaAR A = sPP(J A ) = <0 A A |+<1 A 1 A |. Then 
we can contract these these two morphisms along a single edge as in the following 
picture: 



R A 



La 



This gives us the morphism 10^X0^1 + |1aX1a| which is the identity morphism on 
A. To prove that 2? has daggers, we need to prove a few other things first. □ 

Definition 5.2. The anti-transpose of a matrix N, denoted by N, is N flipped 
across the non-standard diagonal. 

Lemma 5.3. Pf (TV) = Pf(AT). 

Proof. Let ./V = {r/ij} be an n x n matrix. If n is odd, the above is trivial, so let n 
be even. Now let & be the set of partitions of [n] into pairs, (ik,jk), *fc < 3 k- If 
7T e & we can define the sign of n, sgn(7r). This is done by considering the set [n] 
as a sequence of nodes laid out horizontally and labeled 1, . . . , n from left to right. 
Then if two nodes are paired in n, connect them with an edge. Then sgn(7r) is 
(— l) k where k is the number of places where lines cross. Now we can define Pf(iV) 
as follows: _ 

Pf(JV) = ^ sgn(Tr) 

(**,jfc)eir 

Now let rilj = f]n-j+i,n-i+i be the entries of N and suppose 7r e Then the 
mapping & —> & : ir >-» tt' given by (ik,jk) !-» (n — Jfe + 1, n — ife + 1) is a bijective 
involution. Note that n' is the matching formed from tt by relabeling the nodes as 
n, . . . , 1 from left to right. This preserves the number of crossings of edges so that 
sgn(7r') = sgn(7r). Thus we get 

Pf(JV) = J] sgn( 7 r) V' ikjk = 

(ik.jk) 

J] SS 11 ! 71 "') EI Vn-j k + l,n-i k + l = Pf(A f ). 

Tr'eJ? (n-j k + l,n-i k + l) 

□ 

Definition 5.4. If 7 is a bitstring. Let I be the bitstring reflected across a vertical 
axis. Then |7) and (I\ are the tensors where the / is formed from considering I 
as a bitstring. If I c N, I is formed by considering I as bitstring representing a 
characteristic function. Then I is a characteristic function defining another subset 
of N. 

Corollary 5.5. Let N be a skew symmetric matrix with labels M. Let N also have 
labels M. sPf(TV) = £/ QM Pf(JV»|7> 
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Proof. Let J c M. Note that TVj = N f . Then Pf(iV/) = Pf(iVj). This gives the 
result. 

□ 

Proposition 5.6. For any skew- symmetric matrix M, 

2>f(M,)</| 

E pf ( M /)i 7 > 

are morphisms of 2? . This implies that 2? is a dagger monoidal category. 

Proof. Let M have labels A = {A x , . . . , A n ). Then M will have labels A = 
{A n , . . . , Ai). Let Ra be defined as: 



R A = sPf 



where / is the identity matrix reflected over a vertical axis. Then consider the 
following morphism in 



A 


A 





I 


-F 








A n 




Ra 


> 

> 


sPf(M) 


A l 


A 1 




•> — 






! 





This diagram represents the morphism 

S pf ( M ^)i J > 



lei 



S < J i< 7 i 

i^{Aa} 



S Pf(M/)</| 

ISA 

2/Pf(Mj)|7) can be formed similarly by instead using sPf v (M) and sP^i?^). 
Now since every generating morphism has a dagger, the entire category has a dagger 
and it is the usual vector space dagger. □ 

Note that there are two primary types of morphisms in ZP, namely those of the 
form sPf(M) and those of form sPf v (M). The diagrams will then form bipartite 
graphs. 

Suppose we are given a Pfaffian circuit P. Let Si be the gates of the form sPf(M) 
and Qi be the gates of the form sPf v (M). We define S = and & likewise. 

© is the direct sum with the row and columns reordered by a particular ordering. 
The ordering is found by drawing a planar curve through the Pfaffian circuit such 
that every edge is intersected by the curve once and exactly once. Since a Pfaffian 
circuit is planar and bipartite, such a curve always exists and the curve chosen does 
not matter. The edges are then labeled based on when the curve intersects them. 
This is ordering used to define 0. G> is defined to be {(— iy + i +1 6 i3 }. 
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Theorem 5.7. The value of a Pfaffian circuit T is given by Pf(3 + 6)|15j 

This tells us that Pfaffian circuits can be computed in polynomial time. Now we 
want to investigate if there is a functor that transforms determinantal circuits into 
Pfaffian circuits. We want the trace to be preserved so that the resulting Pfaffian 
circuit still solves the same problem as the original determinantal circuit. We also 
want the functor to be faithful. 

Given a morphism in 2?, we show how to construct it in Now it is well know 
that for an n x n matrix M, 



Pf 







M 




det(M) 



We want a way to embed a matrix M into skew-symmetric matrix such that the 
Pfaffian corresponds to the determinant. Let us define M as the matrix M reflected 
across a vertical axis. We define S(M) to be 

Ml 
-M T 



S{M) 



Proposition 5.8. For an n x n matrix M, 





Pf(5(M)) = Pf 



-M 1 



M 




det(M) 



Proof. In general, M can be made from M with [^J column swaps. So if n = 0, 1 
modulo 4, [^J is an even number and so det(M) = det(M). Now if n is congruent to 
or 1 modulo 4, then Pf(5(M)) = (-l)"^- 1 )/ 2 det(M) = det(M) = det(M). If n 
is congruent to 2 or 3 modulo 4, then [^J is an odd number so det(M) = — det(M) 
and Pf(5(M)) = (_i)»(«-i)/2 det(M) = -det(M) = det(M). 

□ 

Theorem 5.9. Every morphism inT> is a morphism in . Thus there is a trace- 
preserving faithful strict monoidal functor from T> — > £P given by inclusion. 

Proof. First suppose that M is an n x n matrix. The labels of S(M) = R u C 
where R is the row labels of M and C are the column labels of M. Now let K be 
a subset of the labels. Then let I = K n R and J = K n C . Then we get 

m z>j 1 





Pf(5(M) K ) = Pf 



-ML 



det(M/,j) 



So we get that 

sPf(S(M)) = ^ det(M 7)J |/>|J> 

sPP(S(M)) ^ det(M r>J )<J|<J| 

Let us consider the identity morphism on A n (g) 
matrix 





A n 


A n -! ■ 


• A, 


A n 


( 1 





' °\ 


A n -x 





1 


■ 


Ax 


1 o 





• 1 / 



<g) Ai in C. It is given by the 
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Suppose we have an n x n matrix M : B\ (g) • • • (g) B n — » A\® ■ ■ ■ ® A n . Then we 
define M* = sPf(5(M)) and i?g>Ai = sPf v (S'(/g )j 4 i )). Let us consider the morphism 
in given by 

Sr 



! 






A n 








' ! 




A l 






• ! 






S 





For J c {Si, . . . , B n }, J, J' c {A„, . . . , Ai}; and J' c [A\, . . . A„}, we can repre- 
sent this tensor as 

/ 2det(M /;J )|7>|J>V2<J'|<J'j 

^]det(M 7i j)|/><J| = sDet(M) 

So for any square matrix M, sDet(M) is a morphism in Now not every mor- 
phism in C is a square matrix. However, if we have an n x m matrix M, we can 
make it square. If n < to, then let M' = M ® Z m - n where Z m _ n is the (to — n) x 
matrix. If m < n, then let M' = M © Z' n _ m where is the Ox (n — to) matrix. 

What this amounts to is either adding rows or columns of zeros as needed. 

Now note that sPf ([0]) = |0>. <0| is also a morphism in Consider sPf v {K) = 
(0a0 b \ + <U1s| where 

A B 

K= A (° 1 
B \-l 

Contracting this with the morphism |0b), we get (0a|- 

Let M be an arbitrary n x to matrix. Then let us consider S(M') where M' is 
defined as above. Suppose n < to. Then 

sPf(5(M')) = 2 det(M /;J )|70 n+ i • • • m ><J| 
i, j 

Consider the following diagram in 



1 — > 

<0 " Z^H sPf(S(M')) 
<0|^ 



1 



TO 



The morphism this represents will obviously come out to be sDet(M). If n > to, 
then copies of |0> are added to the extra output wires of sPf(5(M')). Thus we have 
finished the proof of theorem. Every morphism of T> is in fact a morphism in & . 
Furthermore, the reinterpretation of a determinantal circuit as a Pfaffian circuit 
can obviously be done in polynomial time. □ 
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Despite this fact, detcrminantal circuits still have some advantages. If a Pfaffian 
circuit can be represented as a detcrminantal circuit, its evaluation will be more 
efficient. Also, given the non-intuitive nature of the inclusion, viewing a circuit 
as determinantal may prove more useful than viewing it as a Pfaffian circuit. In 
the applications above, this is certainly true. Finally, the practical complexity of 
implementation is lessened by the determinantal approach since there is less need 
to track a complex ordering on the objects. 
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